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PATENT 

Attorney Docket No.: 018563-004620US 
Client Reference No.: AT00069.2 

EFFICIENT DATA REPRESENTATION OF TEETH MODEL 



CROSS-REFERENCES TO RELATED APPLICATIONS 
[0001] This present application is a continuation of U.S. Application No. 10/241,240 

5 (Attorney Docket No. 018563-004610US - AT-00069.1), filed September 10, 2002, which 
was a continuation of U.S. Application No. 09/506,419 (Attorney Docket No. 18563- 
004600 -AT-00069), filed February 17, 2000, the full disclosures of which are incorporated 
herein. 

BACKGROUND OF THE INVENTION 

10 [0002] 1. Field of the Invention . The invention relates generally to the field of 

orthodontics and, more particularly, to computer modeling of teeth. 
[0 003] 2. Description of the Background Art . Tooth positioners for finishing 

orthodontic treatment are described by Kesling in the Am. J. Orthod. Oral. Surg. 31 :297-304 
(1945) and 32:285-293 (1946). The use of silicone positioners for the comprehensive 

1 5 orthodontic realignment of a patient's teeth is described in Warunek et al. (1 989) J. Clin. 

Orthod. 23:694-700. Clear plastic retainers for finishing and maintaining tooth positions are 
commercially available from Raintree Essix, Inc., New Orleans, Louisiana 70125, and Tru- 
Tain Plastics, Rochester, Minnesota 55902. The manufacture of orthodontic positioners is 
described in U.S. Patent Nos. 5,186,623; 5,059,118; 5,055,039; 5,035,613; 4,856,991; 

20 4,798,534; and 4,755,139. 

[0004] Other publications describing the fabrication and use of dental positioners 

include Kleemann and Janssen (1996) J. Clin. Orthodon. 30:673-680; Cureton (1996) J. Clin. 
Orthodon. 30:390-395; Chiappone (1980) J. Clin. Orthodon. 14:121-123; Shilliday (1971) 
Am. J. Orthodontics 59:596-599; Wells (1970) Am. J. Orthodontics 55:351-366; and 

25 Cottingham (1 969) Am. J. Orthodontics 55:23-3 1 . 

[0005] Kuroda et al. (1 996) Am. J. Orthodontics 1 1 0:365-369 describes a method for 

laser scanning a plaster dental cast to produce a digital image of the cast. See also U.S. 
Patent No. 5,605,459. 

[0006] U.S. Patent Nos. 5,533,895; 5,474,448; 5,454,717; 5,447,432; 5,431,562; 

30 5,395,238; 5,368,478; and 5,139,419, assigned to Ormco Corporation, describe methods for 
manipulating digital images of teeth for designing orthodontic appliances. 



[0007] U.S. Patent No. 5,01 1 ,405 describes a method for digitally imaging a tooth 

and determining optimum bracket positioning for orthodontic treatment. Laser scanning of a 
molded tooth to produce a three-dimensional model is described in U.S. Patent No. 
5,338,198. U.S. Patent No. 5,452,219 describes a method for laser scanning a tooth model 
5 and milling a tooth mold. Digital computer manipulation of tooth contours is described in 
U.S. Patent Nos. 5,607,305 and 5,587,912. Computerized digital imaging of the jaw is 
described in U.S. Patent Nos. 5,342,202 and 5,340,309. Other patents of interest include U.S. 
Patent Nos. 5,549,476; 5,382,164; 5,273,429; 4,936,862; 3,860,803; 3,660,900; 5,645,421; 
5,055,039; 4,798,534; 4,856,991; 5,035,613; 5,059,118; 5,186,623; and 4,755,139. 

! o BRIEF SUMMARY OF THE INVENTION 

[0008] A computer-implemented method generates a computer model of one or more 

teeth, by receiving as input a digital data set of meshes representing the teeth; creating a 
parametric representation of the digital data set; and displaying the computer model of the 
teeth using the parametric representation. 
1 5 [0009] Implementations of the method include one or more of the following. The 

parametric representation is a compressed version of the digital data set. The parametric 
representation is stored in a file. The parametric representation is transmitted to a remote 
computer and displayed on the remote computer. The parametric representation is 
transmitted over a network, which can be a wide area network such as the Internet. The 
20 parametric representation can be a curve network which is fitted to the digital data set. The 
fitting includes obtaining a location for each sample point; obtaining a normal for the sample 
point; generating a (u, v) representation for the sample point; and building the curve network 
from the (u, v) representation. The fitting can include obtaining a location for each sample 
point; obtaining a tangent for the sample point; generating a (u, v) representation for the 
25 sample point; and building the curve network from the (u, v) representation. The method also 
includes selecting a set of sample points from the digital data set; generating tangent values 
for the sample points to define a curve; and determining the length of the curve. 
[0010] Advantages of the system include one or more of the following. The system 

reduces the amount of data storage space required for storing and communicating teeth 
30 treatment information. By utilizing space more efficiently, the system reduces the cost of the 
system, improves the responsiveness of the system , and allows additional functionality to be 
implemented. 
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[001 1] The system also allows visualization to be used to communicate treatment 

information in a computer-automated orthodontic treatment plan and appliance. The 
invention generates a realistic model of the patient's teeth without requiring a user to possess 
in-depth knowledge of parameters associated with patient dental data compression. 
5 Additionally, expertise in 3D software and knowledge of computer architecture is no longer 
needed to process and translate the captured medical data into a realistic computer model 
rendering and animation. 

[0012] The invention allows teeth plan treatment to be generated and communicated 

in a simple and efficient manner. It also improves the way a treating clinician performs case 

1 0 presentations by allowing the clinician to express his or her treatment plans more clearly and 
gives a prospective patients an opportunity to visualize the facial changes associated with the 
proposed treatment. The invention allows multidisciplinary work teams to deal easily and 
efficiently with the treatment plan. Another major benefit is the ability to visualize and 
interact with models and processes without the attendant danger, impracticality, or 

1 5 significantly greater expense that would be encountered in the same environment if it were 
physical. Thus, money and time are saved while the quality of the treatment plan is enhanced. 
[0013] The details of one or more embodiments of the invention are set forth in the 

accompa-nying drawings and the description below. Other features, objects, and advantages 
of the invention will be apparent from the description and drawings, and from the claims. 

20 BRIEF DESCRIPTION OF THE DRAWINGS 

[0014] FIG. 1 A illustrates a patient's jaw and provides a general indication of how 

teeth may be moved. 

[0015] FIG IB illustrates a single tooth from FIG. 1 A and defines how tooth 

movement distances are determined. 
25 [001 6] FIG. 1 C illustrates the j aw of FIG. 1 A together with an incremental position 

adjustment appliance. 

[001 7] FIG. 2 is a block diagram illustrating steps for producing a system of 

incremental position adjustment appliances. 

[0018] FIG. 3 is an illustration of a 3D model of teeth using triangular meshes. 

30 [0019] FIG. 4 is a flow chart illustrating a process for communicating teeth related 

information to a treating professional. 

[0020] FIG. 5 is a flow chart illustrating a process for generating a curve network 

representative of a tooth. 
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[0021] FIG. 6A is a flow chart illustrating a first process for generating an MxN 

network representative of a tooth. 

[0022] FIG. 6B is a flow chart illustrating a second process for generating an MxN 

network representative of a tooth. 
5 [0023] FIGS. 6C and 6D are exemplary diagrams illustrating the operation of the 

process of FIG. 6B on a tooth model. 

[0024] FIG. 7 shows an original high resolution teeth model and various compressed 

versions of the teeth model using the processes of FIGS. 5-6. 
[0025] FIG. 8 shows a diagram of a system for fabricating appliances. 

1 0 [0026] FIG. 9 is a block diagram of a data processing system to develop orthodontic 

treatment plans. 

DETAILED DESCRIPTION OF THE INVENTION 
[0027] Systems and methods are provided for modeling teeth in creating a treatment 

plan. The systems and methods compress data associated with the teeth model to reduce the 

1 5 storage size and transmission time to a treating professional. Based on the treatment plan, the 
teeth can be incrementally moved using a plurality of discrete appliances, where each 
appliance successively moves one or more of the patient's teeth by relatively small amounts. 
[0028] Referring now to FIG. 1 A, a representative j aw 1 00 includes sixteen teeth, at 

least some of which are to be moved from an initial tooth arrangement to a final tooth 

20 arrangement. To understand how the teeth may be moved, an arbitrary centerline (CL) is 

drawn through one of the teeth 102. With reference to this centerline (CL), the teeth may be 
moved in the orthogonal directions represented by axes 104, 106, and 108 (where 104 is the 
centerline). The centerline may be rotated about the axis 108 (root angulation) and 104 
(torque) as indicated by arrows 1 10 and 1 12, respectively. Additionally, the tooth may be 

25 rotated about the centerline, as represented by arrow 1 14. Thus, all possible free-form 
motions of the tooth can be performed. 

[0029] Referring now to FIG. IB, the magnitude of any tooth movement is defined in 

terms of the maximum linear translation of any point P on a tooth 102. Each point Pi will 
undergo a cumulative translation as that tooth is moved in any of the orthogonal or rotational 
30 directions defined in FIG. 1 A. That is, while the point will usually follow a non-linear path, 
there will be a linear distance between any point in the tooth when determined at any two 
times during the treatment. Thus, an arbitrary point Pi may in fact undergo a true side-to-side 
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translation as indicated by arrow di, while a second arbitrary point P 2 may travel along an 
arcuate path, resulting in a final translation d 2 . In many situations, the maximum permissible 
movement of a point Pi in any particular tooth is defined as the maximum linear translation of 
that point Pj on the tooth which undergoes the maximum movement for that tooth in any 
5 treatment step. 

[0030] One tool for a incrementally repositioning the teeth in a patient's jaw is a set 

of one or more adjustment appliances. Suitable appliances include any of the known 
positioners, retainers, or other removable appliances which are used for finishing and 
maintaining teeth positions in connection with conventional orthodontic treatment. As 
1 0 described below, a plurality of such appliances can be worn by a patient successively to 

achieve gradual tooth repositioning. A particularly advantageous appliance is the appliance 
100, shown in FIG. 1C, which typically comprises a polymeric shell having a cavity shaped 
to receive and resiliently reposition teeth from one tooth arrangement to another tooth 
arrangement. The polymeric shell typically fits over all teeth present in the upper or lower 
15 jaw. Often, only some of the teeth will be repositioned while others will provide a base or 
anchor region for holding the repositioning appliance in place as it applies the resilient 
repositioning force against the tooth or teeth to be repositioned. In complex cases, however, 
many or most of the teeth will be repositioned at some point during the treatment. In such 
cases, the teeth which are moved can also serve as a base or anchor region for holding the 
20 repositioning appliance. The gums and the palette also serve as an anchor region in some 
cases, thus allowing all or nearly all of the teeth to be repositioned simultaneously. 
[0031] The polymeric appliance 100 of FIG. 1C is preferably formed from a thin 

sheet of a suitable elastomeric polymeric, such as Tru-Tain 0.03 in. thermal forming dental 
material, marketed by Tru-Tain Plastics, Rochester, Minnesota 55902. In many cases, no 
25 wires or other means are provided for holding the appliance in place over the teeth. In some 
cases, however, it is necessary to provide individual attachments on the teeth with 
corresponding receptacles or apertures in the appliance 100 so that the appliance can apply 
forces that would not be possible or would be difficult to apply in the absence of such 
attachments. 

30 [0032] FIG. 2 shows a process 200 for producing the incremental position adjustment 

appliances for subsequent use by a patient to reposition the patient's teeth. As a first step, an 
initial digital data set (IDDS) representing an initial tooth arrangement is obtained (step 202). 
The IDDS may be obtained in a variety of ways. For example, the patient's teeth may be 
scanned or imaged using X-rays, three dimensional X-rays, computer-aided tomographic 
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images or data sets, or magnetic resonance images, among others. More details on the 
contact or non-contact scanners are in commonly-owned and co-pending Application Serial 
No. 09/169,276, filed October 8, 1998, the content of which is incorporated by reference. 
[0033] A plaster cast of the patient's teeth is obtained by well known techniques, such 

5 as those described in Graber, Orthodontics: Principle and Practice, Second 

Edition, Saunders, Philadelphia, 1969, pp. 401-415. After the tooth casting is obtained, the 
casting is digitally scanned by a scanner, such as a non-contact type laser or destructive 
scanner or a contact-type scanner, to produce the IDDS. The data set produced by the 
scanner may be presented in any of a variety of digital formats to ensure compatibility with 
1 0 the software used to manipulate images represented by the data. In addition to the 3D image 
data gathered by laser scanning or destructive scanning the exposed surfaces of the teeth, a 
user may wish to gather data about hidden features, such as the roots of the patient's teeth and 
the patient's jaw bones. This information is used to build a more complete model of the 
patient's dentition and to show with more accuracy and precision how the teeth will respond 
1 5 to treatment. For example, information about the roots allows modeling of all tooth surfaces, 
instead of just the crowns, which in turn allows simulation of the relationships between the 
crowns and the roots as they move during treatment. Information about the patient's jaws 
and gums also enables a more accurate model of tooth movement during treatment. For 
example, an x-ray of the patient's jaw bones can assist in identifying ankylose teeth, and an 
20 MRI can provide information about the density of the patient's gum tissue. Moreover, 
information about the relationship between the patient's teeth and other cranial features 
allows accurate alignment of the teeth with respect to the rest of the head at each of the 
treatment steps. Data about these hidden features may be gathered from many sources, 
including 2D and 3D x-ray systems, CT scanners, and magnetic resonance imaging (MRI) 
25 systems. Using this data to introduce visually hidden features to the tooth model is described 
in more detail below. 

[0034] The IDDS is manipulated using a computer having a suitable graphical user 

interface (GUI) and software appropriate for viewing and modifying the images. More 
specific aspects of this process will be described in detail below. 
30 [0035] Individual tooth and other components may be segmented or isolated in the 

model to permit their individual repositioning or removal from the digital model. After 
segmenting or isolating the components, the user will often reposition the tooth in the model 
by following a prescription or other written specification provided by the treating 
professional. Alternatively, the user may reposition one or more teeth based on a visual 
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appearance or based on rules and algorithms programmed into the computer. Once the user is 
satisfied, the final teeth arrangement is incorporated into a final digital data set (FDDS) (step 
204). 

[0036] The FDDS is used to generate appliances that move the teeth in a specified 

5 sequence. First, the centers of each tooth model may be aligned using a number of methods. 
One method is a standard arch. Then, the teeth models are rotated until their roots are in the 
proper vertical position. Next, the teeth models are rotated around their vertical axis into the 
proper orientation. The teeth models are then observed from the side, and translated 
vertically into their proper vertical position. Finally, the two arches are placed together, and 
1 0 the teeth models moved slightly to ensure that the upper and lower arches properly mesh 

together. The meshing of the upper and lower arches together is visualized using a collision 
detection process to highlight the contacting points of the teeth. 

[0037] In step 204, final positions for the upper and lower teeth in a masticatory 

system of a patient are determined by generating a computer representation of the masticatory 

1 5 system. An occlusion of the upper and lower teeth is computed from the computer 

representation; and a functional occlusion is computed based on interactions in the computer 
representation of the masticatory system. The occlusion may be determined by generating a 
set of ideal models of the teeth. Each ideal model in the set of ideal models is an abstract 
model of idealized teeth placement which is customized to the patient's teeth, as discussed 

20 below. After applying the ideal model to the computer representation, and the position of the 
teeth is optimized to fit the ideal model. The ideal model may be specified by one or more 
arch forms, or may be specified using various features associated with the teeth. 
[0038] Based on both the IDDS and the FDDS, a plurality of intermediate digital data 

sets (INTDDSs) are defined to correspond to incrementally adjusted appliances (step 206). 

25 Finally, a set of incremental position adjustment appliances are produced based on the 
INTDDs and the FDDS (step 208). 

[0039] Fig. 3 shows one exemplary 3D surface model of the teeth. The surface 

topology of a 3D model of teeth on a jaw can be modeled as a set of polygons of appropriate 
sizes and shapes joined at their edges. The set of polygons defining the 3D object is referred 
30 to as the "model" or "mesh" for the 3D object. In one embodiment, the polygons are 

triangles. In this embodiment, a triangle mesh is a piecewise linear surface with triangular 
faces joined along their edges. 

[0040] Many types of scan data, such as that acquired by an optical scanning system, 

provide a 3D geometric model (e.g., a triangular surface mesh) of the teeth when acquired. 
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Other scanning techniques, such as the destructive scanning technique described above, 
provide data in the form of volume elements ("voxels") that can be converted into a digital 
geometric model of the tooth surfaces. In one implementation, a marching cubes algorithm is 
applied to convert the voxels into a mesh, which can undergo a smoothing operation to 
reduce the jaggedness on the surfaces of the tooth model caused by the marching cubes 
conversion. One smoothing operation moves individual triangle vertices to positions 
representing the averages of connected neighborhood vertices to reduce the angles between 
triangles in the mesh. 

[0041] The triangles in Fig. 3 form a connected graph. In this context, two nodes in a 

graph are connected if there is a sequence of edges that forms a path from one node to the 
other (ignoring the direction of the edges). Thus defined, connectivity is an equivalence 
relation on a graph: if triangle A is connected to triangle B and triangle B is connected to 
triangle C, then triangle A is connected to triangle C. A set of connected nodes is then called 
a patch. A graph is fully connected if it consists of a single patch. The processes discussed 
below keep the triangles connected. 

[0042] Once a 3D model of the tooth surfaces has been constructed, models of the 

patient's individual teeth can be derived. In one approach, individual teeth and other 
components are "cut" using a cutting tool to permit individual repositioning or removal of 
teeth in or from the digital data. After the components are Afreed,= a prescription or other 
written specification provided by a treating professional is followed to reposition the teeth. 
Alternatively, the teeth may be repositioned based on the visual appearance or based on rules 
and algorithms programmed into the computer. Once an acceptable final arrangement has 
been created, the final tooth arrangement is incorporated into a final digital data set (FDDS). 
[0043] Fig. 4 shows an embodiment of a process for communicating the 3D model of 

the teeth to the treating professional. Since realistic models have a large volume of data, the 
storage and transmission of the models can be expensive and time consuming. To reduce 
transmission problems arising from the large size of the 3D model, the system compresses the 
data associated with the model. In the embodiment of Fig. 4, the compression is done by 
modeling the teeth meshes as a curve network (step 220). Next, the curve network is 
transmitted to the treating professional (step 222). The curve network can be sent over a 
local area network (LAN) or a wide area network (WAN) such as the Internet. At the 
receiving end, once the curve network is received, the 3D model is reconstructed from the 
curve network for the treating professional to analyze (step 224). 
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[0044] The curve network generation of step 220 is shown in more detail in Fig. 5 . 

First, the process of Fig. 5 obtains one or more sample points from the meshes of a tooth (step 
230). Next, the process computes tangent values for the sample points (step 232). The 
sample points are converted into u, v-space (step 234). Finally, a grid or curve network is 
5 generated that satisfies the points in u, v-space (step 236). This is done by surface fitting the 
curve network to the original data. 

[0045] In one embodiment of Fig. 5, data points are separated into groups. A U-curve 

with a U-knot vector is interpolated through each separate group. Next, one or more V-curves 
are generated by moving along each of the U-curves for a given knot value. The V-curve is 
10 formed by interpolating through these data points. The V-curves intersect the U-curves at 
points of constant u. A new set of U-curves that are constant in v is obtained by connecting 
points of constant v along the V-curves. A network of curves and its relationship to the 
original data points can then be generated. The only curves are drawn are either curves of 
constant u or constant v, and so the regions marked out by them are rectangular in u, v-space 
1 5 and can be represented at patches which also mark out rectangular domains in u, v-space. A 
patch representation is generated from the network to arrive at a full surface description. 
More details on the surface fitting are discussed in pages 101-1 10 of Alan Watt and Mark 
Watt, Advanced Animation and Rendering Techniques (Addison-Wesley Publishing 
Company, Menlo Park, California). 
20 [0046] Fig. 6A shows a flowchart for a process to create an MxN curve network 

associated with a tooth. First, the process determines M, the number of slices corresponding 
to meridian lines crossing the tooth (step 240). Next, the process calculates tangent values 
for a plurality of planes intersecting with the tooth to define a curve associated with one side 
of a slice (step 242). One or more sample points and sample tangent values for the sample 
25 points are selected (step 244) and the length of the curve is computed (step 246). Next, the 
process of Fig. 6 divides the length of the curve the N and generates points associated with 
the MxN curve network (step 248). 

[0047] Each tooth is modeled as though it were a deformed ball. The curve network is 

a network of M meridian lines and N parallel lines. For each meridian line, a plane is 
30 selected to contain the meridian line. A series of points representing an intersection between 
the plane and the triangle mesh modeling a tooth is determined. The series of points 
represent a curve whose length is then determined. The length of the meridian is divided by 
N. If M is 20 and N is 20, a 20 x 20 grid of points is generated, each point being represented 
in 3-D. In this manner, MxN 3D points can be transmitted in place of the voluminous mesh 
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representation of the teeth. To further reduce the size of the file, the coordinates of each 
point are stored using integer values rather than floating point values. 
[0048] Alternatively, a second process 300 for fitting a curve network to a tooth is 

shown in Fig. 6B, and exemplary diagrams illustrating the operation of the process 300 on a 
5 tooth model are shown in Figs. 6C and 6D. In Fig. 6B, during the process of creating a 

separate mesh for each tooth, a tooth basis is created (step 301). The tooth basis has its origin 
inside the tooth such as in the middle of the tooth. For reference, the z-axis points towards an 
occlusion surface, the y-axis points to a labial tooth surface, and the x-axis points to a mesial 
direction. 

10 [0049] Next, the process 300 selects the number of meridians M. M can range 

between 20-40. Also, the process 300 selects the number of points on each meridian N. N is 
typically ranges between 2/3M to M. A set of M half planes 0, . . ., (M-l) is created: each half 
plane goes through the tooth z-axis, and the plane i forms an angle (2*PI/M)*i with the tooth 
x-axis (step 302). Exemplary half planes are shown in Fig. 6C. 

1 5 [0050] For each half plane, an intersection with the tooth mesh surface is determined 

(step 304. This intersection is a polyline that goes from the lower intersection of z-axis with 
the tooth surface to their upper intersection (points A and B). This line is called a meridian 
line. 

[0051] For each meridian line, a length L is determined and S is determined as a 

20 function of L/(N- 1 ) (step 306). From a starting point, the process 300 proceeds along the 
meridian line and find points pO = A, pi, . . ., p[N-l] = B such that the point p[n+l] is at the 
distance S from p[n] along the meridian. Eventually, an M*N list of points in 3D space is 
generated and saved (step 308). The points can be saved as short integer values rather than 
floating point values. 

25 [0052] The data is saved and sent to a treating professional. To reconstruct a surface, 

the grid of M by N points is used to generate a surface that closely resembles a tooth. 
Although 20 is used in one embodiment, other values for M or N can be 30 or 40. 
[0053] To reconstruct a surface from a M*N grid of points, a Hermit spline F is 

specified that maps a rectangle [0, M]*[0, N] in (u, v) parameter plane into 3D space. The 

30 Hermit spline F is a smooth continued function on the [0, M]*[0, N] rectangle. When 
restricted to any [m, m+l]*[n, n+1] unit square, the spline is a bicubic 3D function. The 
restriction of F onto [m, m+l]*[n, n+1] square is denoted by f (omitting the indices (m, n) 
that should signify that functions fare different on different squares). 
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[0054] To determine a Hermit spline, 3D values F, two derivatives dF/du, dF/dv, and 

the second mixed derivative d2F/dudv in each nod point (m, n) are specified. Thus for each 
bicubic function f on [m } m+l]*[n, n+1] square, 3D values f, two derivatives df/du, df/dv, 
and the second mixed derivative d2f/dudv are generated for the four square corners (m, n), 
5 (m+1, n), (m, n+1), (m+1 , n+1). A unique 3D bicubic function is completely determined by 
those values. The bicubic functions on the next squares are smoothly connected because the 
values and derivatives are the same on the common side. 

[0055] The 3D values F, two derivatives dF/du, dF/dv, and the second mixed 

derivative d2F/dudv in the nod points (m, n) are specified as follows: 
10 F(m, n) = (m, n)-th grid point P(m, n); 

dF/du(m, n) = 0.5*( P(m+1, n) - P(m-1, n) ); 
dF/dv(m, n) = 0.5*( P(m, n+1) - P(m, n-1) ); 
d2F/dudv(m, n) = 0. 

[0056] After the Hermit spline on the rectangle [0, M]*[0, N] is determined, its values 

15 in the more dense grid of points K*L are computed to derive a K*L grid of 3D points. Each 
of the four points (k, 1), (k+1, k), (k, 1+1), (k+1, 1+1) are associated with two triangles: {(k, 1), 
(k+1, k), (k, 1+1)} and {(k, 1+1), (k+1, 1+1), (k, 1)}. In this way, 2*K*L triangles are 
generated that are then rendered as a reconstructed tooth. 

[0057] Fig. 7 shows an original high resolution model 250 of a group of lower teeth 

20 and various compressed models of the high resolution model of the teeth. The original high 
resolution model requires 2.373 megabytes of data. In comparison, a 20x20 model 252 
requires 93 kilobytes, a 30x30 model 254 requires 135 kilobytes, and a 40x40 model requires 
193 kilobytes. The reduction in the size of the file reduces the storage as well as transmission 
requirements while maintaining good quality. At least one orthodontist had noted that the 
25 30x30 model is, for his application, as good as the original high resolution model 250 in 
developing a treatment plan. 

[0058] Although the above parametric techniques can be used to compress the storage 

size of a mesh model of a teeth, other techniques can be used. In another embodiment, 
decimation operations can be applied to the mesh of teeth to eliminate data points, which 
30 improves processing speed. After the smoothing and decimation operation have been 

performed, an error value is calculated based on the differences between the resulting mesh 
and the original mesh or the original data, and the error is compared to an acceptable 
threshold value. The smoothing and decimation operations are applied to the mesh once 



again if the error does not exceed the acceptable value. The last set of mesh data that satisfies 
the threshold is stored as the tooth model. 

[0059] Additionally, the mesh model can be simplified by removing unwanted or 

unnecessary sections of the model to increase data processing speed and enhance the visual 
display. Unnecessary sections include those not needed for creation of the tooth 
repositioning appliance. The removal of these unwanted sections reduces the complexity and 
size of the digital data set, thus accelerating manipulations of the data set and other 
operations. After the user instructs the software to erase unwanted sections, all triangles 
within the box set by the user are removed and the border triangles are modified to leave a 
smooth, linear border. The software deletes all of the triangles within the box and clips all 
triangles which cross the border of the box. This requires generating new vertices on the 
border of the box. The holes created in the model at the faces of the box are retriangulated 
and closed using the newly created vertices. 

[0060] In alternative embodiments, the computer automatically simplifies the digital 

model by performing the user-oriented functions described above. The computer applies a 
knowledge of orthodontic relevance to determine which portions of the digital model are 
unnecessary for image manipulation. 

[0061] In another embodiment, data associated with the 3D model includes features 

can be compressed based on observed characteristics of the tooth surfaces. For example, the 
system can use the curvature of a particular molar between the tips of the cusps and the 
gumline to predict the shape of the roots for that molar. Models of typical root and crown 
shapes also can be used to provide predict the shape of the tooth. The information used to 
arrive at the predicted shapes can be used in lieu of the original data, thus effecting data 
compression. 

[0062] Other feature detection techniques use databases of known cases or statistical 

information against which a particular 3D image is matched using conventional image pattern 
matching and data fitting techniques. One such technique, known as "Maximum a posteriori" 
(MAP), uses prior images to model pixel values corresponding to distinct object types 
(classes) as independent random variables with normal (Gaussian) distributions whose 
parameters (mean and variance) are selected empirically. For each class, a histogram profile 
is created based on a Gaussian distribution with the specified mean and variance. The prior 
images supply for each pixel and each class the probability that the pixel belongs to the class, 
a measure which reflects the relative frequency of each class. Applying Bayes' Rule to each 
class, the pixel values in the input image are scaled according to the prior probabilities, then 
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by the distribution function. The result is a posterior probability that each pixel belongs to 
each class. The Maximum a posteriori (MAP) approach then selects for each pixel the class 
with the highest posterior probability. 

[0063] Once downloaded to a remote workstation, a treating clinician is able to view 

a patient's treatment plan and alter or comment on the plan. The client viewer application is 
implemented in a computer program installed locally on a client computer at the clinician's 
site. In one implementation, the viewer program downloads a data file from a remote host, 
such as a file transfer protocol (FTP) server maintained by the treatment plan designer, which 
can be accessed either through direct connection or through a computer network, such as the 
World Wide Web. The viewer program uses the downloaded file to present the treatment 
plan graphically to the clinician. The viewer program also can be used by the treatment plan 
designer at the host site to view images of a patient's teeth. 

[0064] . 3-D images of various orthodontic views can then be rendered after the 
compressed data of the teeth has been received. In this process, an origin point, or "look 
from" point associated with a camera view is generated. Next, a "look at" point or a focus 
point associated with the camera view is determined. In this system, the line from 
LookFromPoint to LookAtPoint defines the direction the camera is shooting at. Additionally, 
a camera Z vector, or up vector, is determined. 

[0065] Pseudo code for generating various orthodontic views is shown below. With 

reference to the pseudo code, the code defines a boundingbox of one mold (2 arches) which is 
the smallest cube containing the molds geometry. Other settings associated with the 
bounding box include: 

Z_Axis: point from lower to upper, 

Y_Axis: point from inside mouse to front teeth (incisors) 

X_Axis: point from center to left. 

FieldOfView: is the open angle, it corresponding to lens 

HalfFieldOfView: FieldOfView * 0.5 

MoldCenter: Center of the BoundingBox 

X_Length: BoundingBox X dimension 

Y Length: BoundingBox X dimension 

Z_Length: BoundingBox X dimension 

X_MIN: minimum X value of the BoundingBox i.e. right most surface cube X value. 
X_MAX: maximum X value of the BoundingBox 
YMIN: minimum Y value of the BoundingBox 
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YMAX: maximum Y value of the BoundingBox 
Z_MIN: minimum Z value of the BoundingBox 
Z MAX: maximum Z value of the BoundingBox 



[0066] RIGHT BUCCAL OVERJET 

CameraLookFromPoint: 

X = 0.5 * MoldCenter.X + 0.5 * X_Max + 0.25 * 

MAX(Y_Length, Z_Length) / tan(HalfFieldOfView); 

Y = MoldCenter.Y 

Z = MoldCenter.Z - 0.25 * MAX(Y_Length, Z_Length) / tan( 

HalfFieldOfView); 

CameraLookAtPoint: 

X = MoldCenter.X + 0.25 * X_Length; 

Y = MoldCenter.Y; 
Z = MoldCenter.Z; 

CameraUpVector: ZAxis; 

[0067] ANTERIOR OVERJET 

CameraLookFromPoint: 

X = MoldCenter.X; 

Y = 0.5 * MoldCenter.Y + 0.5 * Y_Max + 0.25 * 

MAX(X Length, Z Length) / tan( HalfFieldOfView); 

Z = MoldCenter.Z - 0.25 * MAX(X_Length, Z_Length) / 

tan(HalfFieldOf View); 

CameraLookAtPoint: 

X = MoldCenter.X; 

Y = MoldCenter.Y + 0.25 * Y_Length; 
Z = MoldCenter.Z; 

CameraUpVector: ZAxis; 

[0068] LEFT BUCCAL OVERJET 

CameraLookFromPoint: 

X = 0.5 * MoldCenter.X + 0.5 * X_Min - 0.25 * 

MAX(Y_Length, Z_Length) / tan( HalfFieldOfView); 

Y = MoldCenter.Y; 

Z = MoldCenter.Z - 0.25 * MAX(Y_Length, Z_Length) / tan( 

HalfFieldOfView); 

CameraLookAtPoint: 

X = MoldCenter.X - 0.25 * X_Length; 

Y = MoldCenter.Y; 
Z = MoldCenter.Z; 

CameraUpVector: ZAxis; 
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[0069] LEFT DISTAL MOLAR 

CameraLookFromPoint: 

X = MoldCenter.X - 0.25 * X_Length; 

Y = Y_Min - 0.25 * MAX(X_Length, Z_Length) / tan( 

HalfFieldOfView); 

Z = MoldCenter.Z; 

CameraLookAtPoint: 

X = MoldCenter.X - 0.25 * X_Length; 

Y = MoldCenter.Y; 
Z = MoldCenter.Z; 

CameraUpVector: ZAxis; 

[0070] LEFT LINGUAL 

CameraLookFromPoint: 

X = MoldCenter.X + 0. 125 * X_Length; 

Y = MoldCenter.Y; 
Z = MoldCenter.Z; 

CameraLookAtPoint: 

X = MoldCenter.X - 0.25 * X_Length; 

Y = MoldCenter.Y; 
Z = MoldCenter.Z; 

CameraUpVector: ZAxis; 

[0071] LINGUAL INCISOR 

CameraLookFromPoint: 

X = MoldCenter.X; 

Y = MoldCenter.Y; 
Z = MoldCenter.Z; 

CameraLookAtPoint: 

X = MoldCenter.X; 

Y = MoldCenter.Y + 0.25 * Y_Length; 
Z = MoldCenter.Z; 

CameraUpVector: ZAxis; 

[0072] RIGHT LINGUAL 

CameraLookFromPoint: 

X = MoldCenter.X + 0.125 * X_Length; 

Y = MoldCenter.Y; 
Z = MoldCenter.Z; 

CameraLookAtPoint: 

X = MoldCenter.X + 0.25 * X_Length; 

Y = MoldCenter.Y; 
Z = MoldCenter.Z 

CameraUpVector: ZAxis; 
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[0073] RIGHT DISTAL MOLAR 



CameraLookFromPoint: 



X = MoldCenter.X + 0.25 * X_Length; 

Y = Y_MIN-0.25 * MAX( X_Length, Z_Length) / 



5 tan(HalfFieldOfView); 



Z = MoldCenter.Z; 



CameraLookAtPoint: 
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CameraUpVector: 



X = MoldCenter.X + 0.25 * X_Length; 
Y = MoldCenter.Y; 
Z = MoldCenter.Z; 
ZAxis; 



[0074] 



Once the intermediate and final data sets have been created and reviewed by 



an orthodontist or suitably trained person, the appliances may be fabricated as illustrated in 

1 5 FIG. 8. Common fabrication methods employ a rapid prototyping device 200 such as a 

stereolithography machine. A particularly suitable rapid prototyping machine is Model SLA- 
250/50 available from 3D System, Valencia, California. The rapid prototyping machine 200 
selectively hardens a liquid or other non-hardened resin into a three-dimensional structure 
which can be separated from the remaining non-hardened resin, washed, and used either 

20 directly as the appliance or indirectly as a mold for producing the appliance. The prototyping 
machine 200 receives the individual digital data sets and produces one structure 
corresponding to each of the desired appliances. Generally, because the rapid prototyping 
machine 200 may utilize a resin having non-optimum mechanical properties and which may 
not be generally acceptable for patient use, the prototyping machine typically is used to 

25 produce molds which are, in effect, positive tooth models of each successive stage of the 
treatment. After the positive models are prepared, a conventional pressure or vacuum 
molding machine is used to produce the appliances from a more suitable material, such as 
0.03 inch thermal forming dental material, available from Tru-Tain Plastics, Rochester, 
Minnesota 55902. Suitable pressure molding equipment is available under the trade name 

30 BIOSTAR from Great Lakes Orthodontics, Ltd., Tonawanda, New York 14150. The 

molding machine 250 produces each of the appliances directly from the positive tooth model 
and the desired material. Suitable vacuum molding machines are available from Raintree 
Essix, Inc. 

[0075] After production, the appliances can be supplied to the treating professional all 

35 at one time. The appliances are marked in some manner, typically by sequential numbering 
directly on the appliances or on tags, pouches, or other items which are affixed to or which 
enclose each appliance, to indicate their order of use. Optionally, written instructions may 
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accompany the system which set forth that the patient is to wear the individual appliances in 
the order marked on the appliances or elsewhere in the packaging. Use of the appliances in 
such a manner will reposition the patient's teeth progressively toward the final tooth 
arrangement. 

5 [0076] Because a patient's teeth may respond differently than originally expected, the 

treating clinician may wish to evaluate the patient's progress during the course of treatment. 
The system can also do this automatically, starting from the newly-measured in-course 
dentition. If the patient's teeth do not progress as planned, the clinician can revise the 
treatment plan as necessary to bring the patient's treatment back on course or to design an 

10 alternative treatment plan. The clinician may provide comments, oral or written, for use in 
revising the treatment plan. The clinician also can form another set of plaster castings of the 
patient's teeth for digital imaging and manipulation. The clinician may wish to limit initial 
aligner production to only a few aligners, delaying production on subsequent aligners until 
the patient's progress has been evaluated. 

1 5 [0077] FIG. 9 is a simplified block diagram of a data processing system 300 that may 

be used to develop orthodontic treatment plans. The data processing system 300 typically 
includes at least one processor 302 which communicates with a number of peripheral devices 
via bus subsystem 304. These peripheral devices typically include a storage subsystem 306 
(memory subsystem 308 and file storage subsystem 314), a set of user interface input and 

20 output devices 318, and an interface to outside networks 316, including the public switched 
telephone network. This interface is shown schematically as "Modems and Network 
Interface" block 316, and is coupled to corresponding interface devices in other data 
processing systems via communication network interface 324. Data processing system 300 
could be a terminal or a low-end personal computer or a high-end personal computer, 

25 workstation or mainframe. 

[0078] The user interface input devices typically include a keyboard and may further 

include a pointing device and a scanner. The pointing device may be an indirect pointing 
device such as a mouse, trackball, touchpad, or graphics tablet, or a direct pointing device 
such as a touchscreen incorporated into the display, or a three dimensional pointing device, 

30 such as the gyroscopic pointing device described in U.S. Patent 5,440,326, other types of user 
interface input devices, such as voice recognition systems, can also be used. 
[0079] User interface output devices typically include a printer and a display 

subsystem, which includes a display controller and a display device coupled to the controller. 
The display device may be a cathode ray tube (CRT), a flat-panel device such as a liquid 
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crystal display (LCD), or a projection device. The display subsystem may also provide non- 
visual display such as audio output. 

[0080] Storage subsystem 306 maintains the basic required programming and data 

constructs. The program modules discussed above are typically stored in storage subsystem 
5 306. Storage subsystem 306 typically comprises memory subsystem 308 and file storage 
subsystem 314. 

[0081] Memory subsystem 308 typically includes a number of memories including a 

main random access memory (RAM) 310 for storage of instructions and data during program 
execution and a read only memory (ROM) 312 in which fixed instructions are stored. In the 
10 case of Macintosh-compatible personal computers the ROM would include portions of the 
operating system; in the case of IBM-compatible personal computers, this would include the 
BIOS (basic input/output system). 

[0082] File storage subsystem 314 provides persistent (non- volatile) storage for 

program and data files, and typically includes at least one hard disk drive and at least one 
15 floppy disk drive (with associated removable media). There may also be other devices such 

as a CD-ROM drive and optical drives (all with their associated removable media). 

Additionally, the system may include drives of the type with removable media cartridges. 

The removable media cartridges may, for example be hard disk cartridges, such as those 

marketed by Syquest and others, and flexible disk cartridges, such as those marketed by 
20 Iomega. One or more of the drives may be located at a remote location, such as in a server on 

a local area network or at a site on the Internet's World Wide Web. 

[0083] In this context, the term "bus subsystem" is used generically so as to include 

any mechanism for letting the various components and subsystems communicate with each 
other as intended. With the exception of the input devices and the display, the other 
25 components need not be at the same physical location. Thus, for example, portions of the file 
storage system could be connected via various local-area or wide-area network media, 
including telephone lines. Similarly, the input devices and display need not be at the same 
location as the processor, although it is anticipated that personal computers and workstations 
typically will be used. 

30 [0084] Bus subsystem 304 is shown schematically as a single bus, but a typical 

system has a number of buses such as a local bus and one or more expansion buses (e.g., 
ADB, SCSI, ISA, EISA, MCA, NuBus, or PCI), as well as serial and parallel ports. Network 
connections are usually established through a device such as a network adapter on one of 
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these expansion buses or a modem on a serial port. The client computer may be a desktop 
system or a portable system. 

[0085] Scanner 320 is responsible for scanning casts of the patient's teeth obtained 

either from the patient or from an orthodontist and providing the scanned digital data set 

5 information to data processing system 300 for further processing. In a distributed 

environment, scanner 320 may be located at a remote location and communicate scanned 
digital data set information to data processing system 300 via network interface 324. 
[0086] Fabrication machine 322 fabricates dental appliances based on intermediate 

and final data set information received from data processing system 300. In a distributed 

1 0 environment, fabrication machine 322 may be located at a remote location and receive data 
set information from data processing system 300 via network interface 324. 
[0087] The invention has been described in terms of particular embodiments. Other 

embodiments are within the scope of the following claims. For example, the three- 
dimensional scanning techniques described above may be used to analyze material 

1 5 characteristics, such as shrinkage and expansion, of the materials that form the tooth castings 
and the aligners. Also, the 3D tooth models and the graphical interface described above may 
be used to assist clinicians that treat patients with conventional braces or other conventional 
orthodontic appliances, in which case the constraints applied to tooth movement would be 
modified accordingly. Moreover, the tooth models may be posted on a hypertext transfer 

20 protocol (http) web site for limited access by the corresponding patients and treating 
clinicians. 
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